On the definition of temperature using 
time-averages. 



A. Carati* 
February 2, 2008 



ABSTRACT 

This paper is a natural continuation of a previous one by the author, which 
was concerned with the foundations of statistical thermodynamics far from 
equilibrium. One of the problems left open in that paper was the correct 
definition of temperature. In the literature, temperature is in general de- 
fined through the mean kinetic energy of the particles of a given system. 
In this paper, instead, temperature is defined a la Caratheodory, the sys- 
tem being coupled to a heat bath, and temperature being singled out as 
the "right" integrating factor of the exchanged heat. As a byproduct, the 
"right" expression for the entropy is also obtained. In particular, in the case 
of a g-distributions the entropy turns out to be that of Tsallis, which we 
however show to be additive, at variance with what is usually maintained. 
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1 Introduction 



The problem of defining temperature in non-equilibrium situations is a 
quite delicate one (see for example ref. (Q)). For systems constituted of 
particles, it is usually assumed that temperature should be defined, up to 
a constant factor, as the mean kinetic energy. On the other hand, one 
should take into account that the notion of temperature originates from 
thermodynamics, where the notion of "the particles of the system" has no 
sense at all. In thermodynamics (see refs. (Q), (0)) temperature is defined 
using both the second and the zeroth principle. On the one hand, the 
second principle (e.g. in the Caratheodory formulation) insures that there 
exists an integrating factor of the exchanged heat. On the other hand, since 
actually there exist infinitely many such integrating factors, temperature is 
singled out among them by the zeroth principle, namely by the requirement 
that if the system is put in thermal contact with another one, at equilibrium 
the integrating factors of the two systems and that of the compound one 
should have the same value. 

Thus, when the Gibbs distribution exp(— j3H)/Z((3) is used in statisti- 
cal mechanics, it is first checked that is an integrating factor of the 
exchanged heat, but its identification with temperature requires some more 
work. One has to put the system in contact with another one and to as- 
sume that (at equilibrium) the total system too is described by a Gibbs 
distribution having as Hamiltonian the sum of the Hamiltonians of the two 
components. From this, with some further considerations, one then shows 
that the temperature coincides with f3~ 1 (see for example (Q), chapter 
three.). 

The aim of the present paper is to implement the analogous proce- 
dure when the averaging is performed through time-averages rather than 
through phase-averages with Gibbs' measure. 

For what concerns the time-averages, we follow paper (jjj), the general 
set-up and the main results of which are recalled in Section 2. In Section 3 
it will be recalled how temperature is introduced in the standard thermo- 
dynamic way; moreover the definition of "thermal equilibrium" between 
two systems in terms of time-averages will be introduced, and it will be 
shown that temperature actually exists. Section 4 is devoted to the study 
of the thermal contact between a generic system and a thermal bath (a 
notion that will also be defined), while the thermodynamics in the case 
of the Tsallis ^-distribution will be developed in Section 5. Section 6 is 
devoted to clarifying another point that was left open in paper (0) , namely 
the identification of the quantity a introduced there, which is shown to 
be strictly related to free energy. This result is also used for discussing 
the identification of temperature in the case of the Tsallis ^-distribution. 
Finally, the conclusions follow in Section 7. 
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2 Time— averages 



For a system with phase space .M, suppose a sequence {x n }, x n S M, is 
given, depending parametrically on its first element xo. As a particular 
case, one can think of the orbit generated by the iteration of a map, for 
example the time-Ai map induced by the flow of an autonomous Hamil- 
tonian system. Suppose we are interested in computing time-averages of a 
dynamical variable A(x) (a real function on M.) 

1 N 

i(x )= f - for iV»l, 

n=l 

the number N playing the role of the "final" time, thought of as a fixed 
parameter. One can divide the space M. into a large number K of disjoint 
cells Zj (such that Ai = LiZj), and one has then 

3=1 

where Aj is the value of A at a point x € Zj, and rij is the number of 
times the sequence {x n } visits Zj. It is clear that rij depends on xq so 
that, if a certain probability distribution is assigned for the initial data 
xo, correspondingly rij turns out to be a random variable with a certain 
distribution function, which will depend both on the dynamics (i.e. the 
map) and on the distribution of the initial data. So one can speak in 
general of the "a priori probability P that the cell Zj will be visited a 
number of times rij = n" : 1 

P{nj =n) = fj (n) . (1) 

For the sake of simplicity of the exposition, in paper Q) the following 
hypothesis was introduced: 

Hypotesis 1 The quantities rij are independent random variables, condi- 
tioned by'Y^nj = N . 

This however is not at all necessary, and the computations could have been 
performed without it, as will be shown below. 

From the fact that the occupation numbers rij are random variables, 
there follows that also the time-average A(xq) itself is a random variable, 

1 Notice that in |fj) reference was instead made to the corresponding cumulative dis- 
tribution function Fj(n) = P (rij < n). 
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so that it is meaningful to consider its expectation. Denoting by < • > 
expectation with respect to the a priori distribution, one has then 



N ■ 

3=1 

Now, in statistical thermodynamics one does not deal directly with the a 
priori probability, because it is generally assumed that the time-average 
of a certain macroscopic quantity, typically the energy of the system, has 
a given value, which should play the role of an independent variable. So 
we consider the energy of the system, which we denote by e, and its time- 
average £ = Y2j £ j n j /N , and we impose on the numbers ■ , riK the 
condition 

K 



v ^ e 3 n 3 = U = COnst 
3=1 



The problem of computing the a posteriori expectation < A >u of A given 
U, is solved, in the thermodynamic limit (see reference (5)), by the relation 



<^>^-^^(l + a ) ■ (2) 

3 



where prime denotes derivative, and the function Xj( z ) 1S defined through 
the Laplace transform of the probability distribution function Q by 



exp( Xi (*)) d = X>— /» , (3) 



dcf 

while the parameters 9 and a are determined by the equations 

n = -Ei^(^ + «). 

In terms of the quantities 



(4) 



v i = - x AN +a ) ' (5) 
relations J5J and Q take the form 

3 3 3 

and this shows that Vj can be interpreted as the mean occupation number 
of cell Zj . 
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In particular if the process of occupation of any cell is a Poisson one, 
i.e. if the successive visits of a given cell are independent events, then one 
finds 

Xj(z) = pexp(-z) -p , (7) 
with a parameter p > 0. In such a case one easily shows that the system 
follows a Gibbs statistics. In fact the mean occupation numbers are easily 
calculated from (JjJJ, and turn out to be given by 

e -6e k /N 

= N ^W ' (8) 

where e" = Z{9) = J2k e^ 6ek ^ N is the usual partition function, so that 
relations JSJl become the usual canonical ones. 

Let us now recall how thermodynamics was formulated in reference (0) 
in terms of time-averages. First of all, following Gibbs and Khinchin (see 
refs. (0),(0)), one defines the exchanged heat SQ in terms of the work per- 
formed by the system. Indeed if one defines the macroscopic work <5W 
performed by the system as 8W =< d K H >u dn, i.e. as the expectation of 
the microscopic work performed when a parameter n entering the Hamil- 
tonian is changed, then, using the first principle one defines the exchanged 

heat as SQ = f dU — SW. It is easily shown (see ref. (0), but also below) 
that one has 

i 

SjdPj . (9) 



N ^ 

3 



This expression shows, recalling 101, that it is convenient to introduce as 
an independent variable, in place of z, the quantity Vj — — Xj(z), and this 
naturally leads to introducing in place of \j its Legendre transform hj, 
defined as usual by 

hj(vj) = ("jZ + Xj(z))\ v . = _ x , (z) ■ 

Notice that, while Vj has the meaning of a mean occupation number (con- 
ditioned on [/), the quantity Uj just plays the role of a parameter, in the 
same sense as z does in (jSJ). In particular, the quantities Vj do not need 
satisfy any condition related to normalization, or the fixing of an energy 
value. Now, from the Legendre duality, one has 

*i = -X'j (jf + «) + <* = h J (^) . 

so that, expressing Sj in terms of hj and using ^ dvj — 0, relation © 
takes the form 
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This shows that the exchanged heat always admits an integrating factor. 

The problem left open in reference (|5|) is that there actually exist in- 
finitely many integrating factors, so that a further requirement is needed 
in order to single out which one should be identified with the inverse ab- 
solute temperature. The aim of this paper is to show that, under suit- 
able hypotheses, 9/N indeed is the inverse temperature, and consequently 
h = J^j hj(^j) is the thermodynamic entropy. 

3 Thermal equilibrium. 

As a general fact it is well known that, if the exchanged heat SQ admits an 
integrating factor 9 (i.e. 9SQ = da for a certain a), then any function of the 
form 9F{a) will be an integrating factor too. The (inverse) temperature p 
is singled out by the requirement that, if two systems are put in thermal 
contact, at equilibrium the values of the integrating factors are the same. 
More precisely, consider systems A and B, with their integrating factors 
9\, 9i and related functions o\ and &i. Put them in thermal contact to 
form system C, with integrating factor 9\2 and function a\2- Then one can 
show (as recalled in the Appendix) that for each of the systems there exists 
essentially a unique integrating factor, characterized by the property that 
the values of the three functions P2, and (3\2 are equal when A and 
B are in mutual equilibrium. By definition, the inverse of this integrating 
factor is the (absolute) temperature of the system, and the corresponding 
functions Si, S2, S12 are the entropies. From the equality of the integrating 
factors there follows immediately that one also has 

dSi2 — dSi + dS 2 , 

i.e., as one usually says, entropy is additive. We want to implement now 
this thermodynamic approach in order to identify the temperature of our 
system. So we have to couple our system to another one, define the notion of 
mutual equilibrium and show that three corresponding integrating factors 
can be found such that at equilibrium their values are equal. 

Consider two systems with phase spaces A4i and M2, let {Zj} and {Z^} 
be the corresponding partitions into cells, and e*. and ejj, the correspond- 
ing values of the energy. When the systems are isolated, they are supposed 
to be described as in Section 2. When put in thermal contact, they form 

a compound system with phase space M12 = Mi x M2, and {Zjk} = f 
{Zj x Zk} is a partition with corresponding energies + + Ae^*; 
here, the term Xejjf corresponds to the interaction due to the thermal con- 
tact, which should be thought of as small (one may assume A < 1). We 
denote by rijk the numbers of times an orbit {(xn\x^)} in A4i2 visits 
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the cell Zjk, and correspondingly we denote by rij and rik the number of 
times the corresponding projections {x„ ^}, {x„ ^} visit the cell Zj and Zk 
respectively. One obviously has the relation 2 

k 

n k = ^2 n jk ■ 

3 

Concerning the a priori probability distribution for the occupation num- 
bers, this will in general depend not only on the distribution of the initial 
data in the space M.\ x M2, but also on A (i.e. on the dynamics). One 
should however take into account that, if the two systems are at the same 
temperature, nothing happens when they are put in contact, i.e. the prob- 
ability distribution for rij and n k will not change, or rather will change 
so little that the changes can be neglected. This in turn implies that the 
probability distribution on the product phase space cannot be given in 
an arbitrary way, but must have some relation to the case A = 0. With 
this motivation in mind we give the following definition of mutual thermal 
equilibrium 

Definition 1 Two systems are said to be in mutual equilibrium if their 
a priori probabilities ff~\n) and fj^\n) do not depend sensibly on A for 
A~0. 

It can be shown (see the next section) that the notion of thermal equilib- 
rium implies that the mean energies U\ and U2 of the two systems cannot 
be given at will, but (having fixed all external parameters) the value of 
the energy of any of the two systems fixes the value of the energy of the 
other one. Thus, in the plane U\, U2 there remains defined an equilibrium 
curve which determines the relation that the energies of two systems have 
to satisfy when they are in mutual equilibrium. 

To show that there exist integrating factors which have the same value 
for the three systems, one needs taking into consideration an apparently 
obvious relation among the exchanged heats, namely 

SQ12 = SQi +SQ 2 , 

where SQi, SQ2 and <5Qi2 are the heats exchanged by system 1, system 2 
and the compound system respectively. This relation is actually far from 

2 For the sake of notational simplicity, here and in the following we write rij in place 
of n^p , and in place of , with the understanding that index j refers always to 
system 1 while index k refers to system 2. The same is understood for other objects 
such as for example Zj, Z^. 
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trivial, because the single terms 8Q\ and 5Q 2 are not a priori the same 
ones as one would have in the absence of a thermal contact between the 
systems. The relation is however true when there is a mutual equilibrium, 
because in such a case one has 

sq» = 5>? J + 4 2 V% = E # 5 E ^ + E 4 2) E ^ 

jfe j k k j 

= E e ? + E e k )dDk = 6 ®i + ^ 2 ' 

where Djj- is the expectation of rijk conditioned by ^ n^fe = N and XX £ j + 

el )fijk — U. In the above expression, in virtue of Definition ^ the quanti- 
ties Pj, &k are essentially the same as in the uncoupled case A = 0; by the 
same token the contribution \e"£ to Ejk was neglected. From the above ex- 
pression for SQ12 and relation (|10|l there follows that there exist functions 
02, and correspondingly eri, 02, 012 f such that 

— do"i2 = -z-dvi + ^~da 2 ■ (11) 

"12 til t) 2 

From this it can be proved that there exist three integrating factors (the 
inverse temperatures) /3i, f3 2 , P12, uniquely defined apart from a multiplica- 
tive constant, which have the same common value for the three systems, 
and correspondingly three functions (the thermodynamic entropies) Si, S2, 
S12 which are additive in the sense that 

dSi2 = dSi + (IS2 ■ 

The proof is standard and is recalled in Appendix A. 

There remains the problem that, for the coupling of two generic systems, 
we are presently unable to find an explicit expression for the entropy of the 
compound system. We are able however to do it when one of the systems is 
a heat bath. In the next section the coupling of a system with a heat bath is 
considered, and it is shown how to compute the probability distribution of 
the occupation numbers njk and the thermodynamic quantities of interest. 

4 A system in contact with a heat bath. 

We want to compute the probability distribution of the occupation numbers 
rijk for the cells Zjk of the compound system M.\ x AI2, when the second 
system is a heat bath, i.e. follows a Gibbs distribution. For the sake of 
simplicity we will limit ourselves to the case in which the probability /j 1 ' (n) 
for the occupation number of the cell Zj of the first subsystem does not 
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actually depend on j, i.e. one has /- \n) = f^(n). Coherently, the 
corresponding Laplace transform J3j) will be denoted by exp(%(z)). From 
Definition n° nc can limit oneself to consider the uncoupled case A = 0. As 
recalled above, in the case of a system described by a Gibbs distribution the 
visiting of the cells are independent events having a common probability p 
to happen. 

Our aim is now to compute P({rijk}), namely the probability of a given 
set {rijk} of occupation numbers. The main difference with respect to the 
case considered in ref. (0), is that now the random variables rijk cannot be 
assumed to be independent, so that now P ({rijk}) is not factorized. One 
can proceed in the following way. For a given set {rijk} let lj = X)fc n i fe 
be the corresponding number of visits of cell Zj in the first system. As 
the visits of the cells of the second system are independent events, the 
probability distribution conditioned by fixing lj will be multinomial, i.e. 
will be giving by 



k 



On the other hand, if f^\lj) is the probability that the cell Zj of the first 
system is visited lj times, as the occupation numbers of the first system 
have been assumed to be independent (Hypothesis one finally has 

= J] lp l f^Hh). (12) 

In computing the conditional expectations, it will be seen in the next pages 
that an essential role is played by the Laplace transform of distribution 
A simple computation shows that 

e -S* WP({ % ,}) = exp [^x(-log(p^e- 2 *))] , (13) 

{rijk} 3 k 

where by X}{n- fc } we mean a sum over all possible sets {rijk}- In fact, (|13() 
follows from the chain of identities: 

e«- e ^^w)= e e n' j - ;...:. K , m 

{tj>0}{L k n jh =h} 3 3 J 2 ' 

x , ljl(pe- z i 1 ) n i 1 ■ ■ ■ (pe~ z i K 2 ) n J K 2 



II E ^ nn\---n jK2 l 



=nE^)(E^ 

j (j>0 k 



ex p [E^( _1 °g^E e Zjfc )) 



3 
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One has now to compute the expectation of the time-average of a 
generic dynamical variable A, conditioned by XX £ j^ + e i 2 ^) n jfc = ^ an( i 
^Tijfe = N, namely the quantity 

< A >u= ^Yl'Jl A ^ k P ^ n Jk)) I Yl p (i n ik}) > 

{«Jfc} i k {n jk } 

where ^ denotes a sum over the possible sequences {njk} constrained 

by J2 n jk = N and XX £ j + £< k^) n jk — U. This can be reduced to the 
computation of the "generating function" 

Z(A, (j,) d ^ Y! ex V (-t,J2 A ^n jk )P({n jk }) , (14) 



through the relation 

1 F) 

(15) 



It turns out that, as in reference (ja), the asymptotic expansion of the 
generating function Z(A, fj,) is very simply computed in the limit of very 
"large" systems (the ones of interest for thermodynamics), by using the 
steepest descent method. This indeed is commonly done in statistical me- 
chanics, following Fowler and Darwin (see ref. (7)). Such an expansion will 
be performed below up to the leading term, neglecting the remainder (an 
explicit expression of which could however be given). As the remainder 
depends both on the form of the energy of the total system, and on the 
function x( z ) (i- e - on the probability distribution) characterizing the first 
system, the validity of the procedure should be checked for any particular 
system. In paper (|5|) it was shown that such a procedure is indeed correct, 
for example, for systems described by the Gibbs measure. So, we suppose 
that our systems too are well described by the leading order term of the 
asymptotic expansion, and we presently show how the expansion is actually 
performed. 
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This goes as follows. The expression (|14f) can be rewritten as 

{n ]k >0} 

cx P(-A*E^ fen J fe ) P ({ n J fe }) 

L 

= lim / / dK\d,Ki exp(— in-JJ — iK%N) 

L^+ocJJ 
-L 

E e M-J2 n jk(M 3 k + ^-(4 1} +4 2) ) +»«2)-p({"jfc} 

{n jfc >0} 
L 



lim // dnidK,2exp(—iniU — in2N) 

L— >+oo ' 



exp 



]T x ( - i og (p^ e M^+i«i(4 1, +4 2) )/w+'«= 



where in the second line the familiar representation of the Dirac delta 
function 5(x) = J dnexp^inx) was used, while in the third line use was 
made of formula for the Laplace transform of the probability P({n,jk})- 
The (double) integral in the last line can be evaluated using the steepest 
descent method, and to leading order one finds 



\ogZ{A, l x) = --U-aN- 



e MA 3fc +0(e< 1) +4 2) )/JV+ce 



Ex(- lo s(^E 

3 

where 6 and a are the solution of the system 

u = -jf E(4 X) + 4 23 )x' f - 



exp 



E^ ( / -log(p^e^+^ 1, +4 2 

jk V k 

exp(M jfc + |( £ «+4 2 ))) 



Efc ex P (Mi* + #( £ i 1} + e 



n 



Now, taking the derivative of log Z(A, fi) and putting /i = 0, after some 
simple algebra one finds 

< a >^-mE a ^ + a + log(pZ2) ) "W~ ' (16) 

jk 

where we have defined Z 2 (9) = f ^2ey^p{—9e^/N), whereas the constants 
9 and a are the solution of the previous system with /j = 0, i.e. are solution 
of 

j re 



7V = -E^'(iv : ^ + a+logpZ2) 



4 2) /iv 



(17) 



, 'JV J °' " Z 2 (9) 

jk 



The formulae and (|17fl solve the problem of computing the conditional 
expectation of a dynamical variable of the compound system, when the 
second one is a heat bath. 

Now, it is very interesting to consider two limit cases: that in which 
the observable A depends only on the variables of the first system (so that 
Ajk — Aj), and that in which A depends only on the variables of the 
second one (Ajk — Ak). In the first case it is meaningful to consider a 
situation in which the energy conditioning is given not on the total energy 
of the compound system, but on the energy U\ of the first one. This essen- 
tially amounts to considering the first system as isolated from the second 
one. Now, if one computes the generating function with the conditioning 
Yl e j n jk/N = U\, one finds 



<A> Ul =-±Y J A ^t%4 ) + a ^) > ( 18 ) 

3 

with 9\ and o,\ solution of 



3 



(19) 



These relations are the same as (|2J and J2J of Section 2, as it should be. 

Instead, if the dynamical variable A depends only on the variables of 
the heat bath, and the conditioning is done only on the energy U 2 of the 
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latter, one finds 

< A~> - 1 V A e M-d2e { ^/N) 

< a > U2 - -^\ Ak — z^m — ' (20) 

where Z 2 is the canonical partition function defined above, while 9 2 is 
defined by 

These are the standard Gibbs relations. 

It is interesting to compare these results with the computation of the 
mean energy of any of the two systems, when the total energy U is fixed. 
One finds (using in the second line the definition of Z 2 ) 



e ~9e { k 2) /N 



N ^—f J N 3 a ^ " Z 2 (9) 

jk 

(22) 

j 

which is the same as (11911 with U\ — U\{9) = f < eW >u- For what concerns 
U 2 one finds instead (using now, in the second line, the second of (CHJ) 

1 ^ (2) ,.0 (1) , „ ,e-< 2) l N 



< & >^ = -4fE CVC^f + a + lo gP Z 2 ) 



N ^— ' N - 1 °^ " Z 2 {9) 

jk 

= -s^+'+'wwE'f'w (23) 

j k 

Y k Z 2 (0) ' 

namely again the same as (EU, with U 2 = U 2 {6) = < e( 2 ) >u. These 
computations show first of all that 

U(9) = U 1 (9) + U 2 (0) , (24) 

but also, as mentioned in Section 3, that the equilibrium energies U± and 
U 2 lie on a curve, i.e. the curve (Ui(0), U 2 {9)) parametrized by 9. 
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5 The thermodynamics. 



We turn now to formulae Ijltjfl and Q17[l. in order to write them in a more 
transparent way. In fact, defining in perfect analogy with JSJ the mean 
occupation numbers Vjk by 

- dcf ,(6 (1) \ exp(-fej 2) /AQ 

such formulas take the form 

< A >u = n 52 A ^ k > u = n 52(4 ] +4 2) )^ . ^ = 52 v * ■ 

jk jk jk 

As in the case of (jSJ, one has now 





Vjk 



77 — ( -logp Vexp(-z r; ) ] 



2,/,=(e/JV)( e( ( , 1) + E [ a ')+c 



One can then introduce the Legendre transform /i^ 12 ^(^n, ■ • ■ ) of the 
function x(- • ■ ) occurring above, by 



jk I' \ I 



eM-zin)J , (26) 
where as usual the dependence of Zjk on i/jk is obtained by solving 

Vjk = d Zjk ^2x (^-logpY^exp(-ziu)j ■ (27) 
Now, the Legendre duality gives 

so that for the exchanged heat 5Q n one finds 

jk jk x J 

This shows that the exchanged heat does indeed have an integrating factor. 

But there is more. Indeed, from the definitions JSJ and (|25|l of the 
quantities Vj and Vjk respectively one checks that 



52 



Vjk 
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so that 

^E^ = E^ = ^ = y^ (1) > 

jk j 

where the last equality comes from (I10II . In the same way one can check 
that 

^E4 2) ^=E4 2) ^ = ^ = ^ (2) - 

jk k 

Here is the standard Gibbs entropy ft/ 2 ) = ^ n k logrifc — logp, while 
Dk is given by ijHJ- 

Finally, we find that 

^dh^ = SQ 12 = SQ 1 + SQ 2 = %dhW + ^dh^ . 

u u u 

This shows that 6/N is indeed the absolute temperature, while the ther- 
modynamic entropies can be identified 3 as 

S 1 (U 1 ,k) d ^ hWfo) , S 2 (U 2 ,k) d ^ f hW{D k ) , S 12 (U,k) d ^ f h^ 12 \v jk ) . 
Moreover the thermodynamic entropies are additive, in the sense that 

dS\ + dS 2 = dS\2 ■ 

6 The case of the q— distribution. 

If system 1 is described by a Gibbs statistics, the function \ 1S an expo- 
nential, and the usual formulae are recovered. In fact in such a case one 
gets 

Ex(-io gP (E^(4 1) +4 2) )+«)) =E^p(->f } +4 2) H • 

j k jk 

From this, the additivity of entropy easily follows using the formula 

h = y^^v jk \og9j k , 

jk 

and the property 

1 _ _ 

3 Notice that in paper 0) the thermodynamic entropy was denoted by S th , while the 
quantity h was denoted by S. 
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which follows from definition In fact, one has 



P jk =pe- a e- e ^/ N e- e ^/ N 



On the other hand, computing a by the second of l|17(l . one finds 



pe a e 



Zi(9)Z 2 (0) ' 



so that for za^ one gets 



Vnh = N — — = ViVh . 

Jk Z X (Q) Z 2 {6) N 

the second equality following from the known expression JSJ which holds 
both for T>j and i>k- 

A comment is now in order before discussing the expression of the en- 
tropy in the case system 1 follows the Tsallis ^-distribution. The relation 
Djk = jjVjVk, which can be rewritten in the more expressive form 

< rijk >u(0)= < rij >u 1 (8)< n k >u 2 (e) , (28) 

holds always true, as one can check from the definition l|25|l of Vjk and 
from the definitions JSJ) and iJSj of v.j and T>u- But the important point 
is that a relation of this kind holds neither for the occupation numbers 
n,jk, Tij, 7ifc, nor for the parameters i/jk, Vj, Vk- In particular, concerning 
the two functions h^ 12 \vjk) and h^^Vj) + h^(vk) (the values of which 
obviously cannot be compared unless some relation is assumed between 
the corresponding arguments), we will find that their values coincide when 
they are computed at equilibrium, i.e. when their arguments satisfy relation 

(23). 

If system 1 follows the Tsallis g-distribution, then one has 

X (z)=Pi(l-(l~q)z) 1/{1 - q) -Pi (29) 
with a constant pi > 0, so that (from EE) one gets 

v jk = Pi (1 + (1 - q) logpZ 3 ) q/{1 - q) C -^- (30) 

with Zj = X)fc e ~ Zjfc - The aim is now to compute the explicit expression 
of the function ft/ 12 ', the entropy of the compound system, and make a 
comparison with the sum + h^ 2 \ 
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First one has to express Zj k as a function of Uj k . To this end, note that 
from l|30|l . taking the logarithm of both sides, one gets 



Zjk 



logz/j& - log ( 1 - (1 - q) logpZj 



9/(1-9) 



log Pi -logZj 



In turn, the expression of the function Zj in terms of Vj k is obtained from 
H30fl by summing over the index fc, which gives 



^2 v i h = Pi C 1 + C 1 - V) 1°EP Z j 



,9/(1-9) 



so that 



log pZj 



V pi / 



Inserting this relation into the expression for Zj k one hnds 



jfc = ~ log I/,-* + log ( ^fc) + J" - 



E l — 9 
__k^jk\ i 

Pi 



\ogp . 



It is now immediate to perform the Legendre transform M 12 ) of the 
function J2jX{-) = J2jPi( l ~ 0- ~ ?) logO^)/?) 17 ^ 9 ' - Pi, and one 



obtains 



h (i2) = _^(- Vjk log v jk + v jk logp) + XI ( 51 ^ fc ) log ( H 



1 Pi9 V 



Pi 



This expression reduces at equilibrium, i.e for i>j k = Vj v k /N, to the simpler 



h (12) = ^2^ Vk \ ogVk + v k \ gp) +p v 



' ^3 ( Pi ) ( Pi 



The first term coincides with the familiar expression of the Boltzmann 
entropy of the heat bath, whereas the second one coincides with the ex- 
pression of the Tsallis entropy for the first system (see refs. (0),(0)), written 
however in terms of the mean occupation numbers, rather than in terms of 
the escort probabilities. It is then apparent that 



dh^ = dhW + d/i< 2 > 
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7 The meaning of the parameter a. 

We now come back to another problem concerning the case of a single 
system, which was left open in paper |Sj), namely the identification of the 
quantity a entering formula (J2J. To this end we remark that, according to 
formula Ipfl. the mean occupation number is given by 




whereas the parameter 8/N was shown to be the temperature of the system. 

It is thus convenient to choose as independent variables (3 = 9/N and N 
in place of U and N. In fact if we define 

/(/?) = -^ExiG^ + o). 
j 

taking the derivative of both sides with respect to j3 and using JfjJ) one finds 
(we denote ^ — dp) 

d f3 (-a-f) = U . 

Now we recall that the free energy F = U — /3 _1 S satisfies the relation 

d (0F) = U , 

so that one has, apart from an additive constant (possibly depending on 
the external parameters entering the Hamiltonian) , the important relation 

a=-pF-f. 

In particular, if the system follows a Gibbs statistics, / is a constant 
and one has 

a = -/3F . 

This coincides with the familiar relation (3F = — logZ, where Z is the 
canonical partition function, because in such a case one has a = log Z (see 
formula JBJ). 

If instead the system follows a Tsallis distribution of index q, then the 
expression of xi z ) is given by l|29() . so that / is given by 

^EM 1 -! 1 - 9)(J3ej+a))^ -pi] 
j 

= Jf I> i 1 - t 1 - q W £ * + ")) i 1 ~ ( 1 - q ^ £ i + «)) ^ - 1 
j 

= !-(!- q)(/3U + a) - 1 = (g - l){flU + a) , 
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where, in the second line, use was made of the second and the third of 1(5} . 
We find in this way 

a=~0[F-(l-q)U] , 

which gives the relation between a and the thermodynamic functions F 
and U. 

At this point, it is worth noting that usually the Tsallis distribution is 
written in the form (see ref. fU)) 

v J = C<j3 q )(l-/3 qBj )T% , 

where the function C((3 q ) is a normalization factor, and (3 q is determined 
by somehow fixing the mean energy. But in terms of our f3, i.e. of inverse 
temperature, one has 

(1 - g)/3 

This shows that j3 q is not the inverse temperature, but a complicated func- 
tion of it, which could be obtained by expressing a as a function of /?. 



8 Conclusions. 

So we have shown, that the parameter /3 = 9/N is the inverse tempera- 
ture, in the sense that it has the same value for every system which is in 
thermal equilibrium with a heat bath at (inverse) temperature (3. Such an 
identification also enables one to find out the thermodynamic entropy S. 
In particular, in the case of the Tsallis g-distribution, the entropy just co- 
incides with his g-entropy (as one could have imagined). A relevant point 
is however that in such a case, at variance with what is usually maintained, 
the entropy turns out to be additive, at least for what concerns its differ- 
ential, i. e. in the sense that one has dSn = dSi + dS^- Whether such a 
relation can be integrated to give S12 = Si + S2 is a non trivial point which 
we are unable to discuss at the moment. In the literature there is a long 
debate about this point; in particular it is often pointed out that, if entropy 
is assumed to be additive (for independent systems), then the Boltzmann- 
Gibbs expression should follow (see ref. ffij)- This is usually based on a 
uniqueness theorem of Khinchin (see ref. I|l(ll ) in the information theory 
framework, in which additivity plays a key role. 

Without entering these very interesting questions, we only want to point 
out that, in our definition 126fl of the entropy, the independent variables are 
the parameters Vj , which are not the probabilities of the occupation of the 
cells in the space phase. So the functional dependence of the entropy on 
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the quantities Vj may be out of the reach of Khinchin's theorem, in which 
the role of the independent variables is played by the probabilities. In any 
case, it is true that one ought to understand in a deeper way the connection 
of the present approach with Kinchin's theorem, and more generally with 
information entropy. 



Appendix: proof of the existence of tempera- 
ture. 

In the following, for the sake of completeness, the familiar deduction of 
the existence of absolute temperature is recalled. The only difference with 
respect to the treatments of most textbooks, is that we make use not of 
the concept of "empirical temperature", but rather of the concept, recalled 
above, of the "equilibrium curve" in the plane {U\, Uz). 

So let us assume there exist 9\, u\ such that for the exchanged heat 
8Q\ of system A one has SQi — 6\dcr\, and analogously for system B and 
the compound system AUB, i.e. 5Q2 — Qidoi , SQ12 = Q\id(J\i- From 
SQ12 = SQi + SQ 2 , one has 

Oi 2 da 12 = Older 1 + 2 da 2 , (31) 



i.e. 



C12 @12 

This shows that 012 depends only on the values of e>i and e>2 , and not on any 
other external parameters entering the energies of system A or system B, 
i.e one has 

C12 = G((7i, a 2) ■ 

So there follows that the ratio of 9i by 612 is equal to the partial derivative 
d ai G, i.e. one has 

9\ = 7; — ^12 , 9 2 = 7; — 6*12 . (32) 
OUi 002 

Use now the condition that systems A and B are in mutual equilibrium. 
This means that, if one fixes all the parameters but the energies Ui, U2 
of the two systems, then they lie on a curve (Ui(p), ^(p) in the plane 
(Ui, 1/2)- For example, in the case of two gases, one can fix the volumes 
and think of changing the internal energies of the gases through an isocoric 
transformation, i.e. by heating or cooling the gases. In our case, instead, 
one can think of fixing the entropy and changing U by adiabatic transfor- 
mations. In other terms, one uses as independent variables the entropies 
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and the internal energies (and, perhaps, some further parameters if the 
former are not sufficient to the complete thermodynamic description of the 
systems) . 

Now, taking the logarithmic derivative of the expressions with 
respect to p, i.e. the variable which parameterizes the equilibrium curve, 
one finds 

dp log 6»i = d p log 6>i2 = dp log 6» 2 . 

In principle, here, <9 p log#i depends only on p and on the other parameters 
of the first system, while d p log 02 depends only on p and the parameters 
of the second one. This implies that actually all the three expressions have 
to be equal to a function depending only on p (as can be seen, for example, 
by the fact that the derivatives with respect to the parameters other then 
p vanish). One has thus 

dp log 6»i = d p log 6i2 = dp log 2 = f{p) , 

which on integration gives 

0! = e F W£i(<7i) ) 2 = e F WEi(<T 2 ), 12 = e F ^£i 2 (<7i, a 2 ) , (33) 

where F{p) is a primitive of f(p) and are integration constants. One 
can wonder why depends only on 0{. This too follows from (|32|l . which 
shows that the ratio 0\/0 2 depends only on ax, cr 2 , and not on any other 
external parameters needed to describe the system. 

The last step is to show that T -1 ^ f exp(F(p)) is the integrating factor 
one is looking for. By construction, T _1 has the same value for all the 
three systems. It is uniquely defined apart from a multiplicative constant 
depending on the choice of the primitive F. To show that it is an integrating 
factor, one has to consider relation l|31|) . which, after simplifying T, reads 

Ei2<icri2 = SidfTi + Y, 2 da 2 ■ 

Thus one can define the entropy of the two systems as dSi = Y,i{a\)da\ 
and dS 2 — Y,2{(?i)d<j2 respectively. Moreover, as one has 

Si 2 dCTi2 = dS\ + dS 2 = d(Si + S 2 ) , 

one finds in the first place that Tii 2 da\ 2 too is a total differential dS*^, and 
furthermore that the entropies are additive in the sense that 

dS\ 2 = dS% + dS2 ■ 

This also shows that T~ l is an integrating factor, because one has 

5Qi = TdSi , SQ 2 = TdS 2 , SQ i2 = TdSi 2 . 
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